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Abstract 

We present a detailed study of the properties of the phase transition in 
the four-dimensional compact U(l) lattice gauge theory supplemented by a 
monopole term, for values of the monopole coupling A such that the transition 
is of second order. By a finite size analysis we show that at A = 0.9 the critical 
exponent is already characteristic of a second-order transition. Moreover, we 
find that this exponent is definitely different from the one of the Gaussian 
case. We further observe that the monopole density becomes approximately 
constant in the second-order region. Finally we reveal the unexpected phe- 
nomenon that the phase transition persists up to very large values of A, where 
the transition moves to (large) negative (3. 
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1. Introduction 



The continuum limit can be recovered from the lattice regularization of a quantum 
field theory by approaching a critical point with suitable scaling properties. While 
the existence of such a critical point can be established through perturbative consid- 
erations in the case of asymptotically free theories, very little is known about critical 
points that can give origin to continuum limits in the non-perturbative domain. The 
£7(1) lattice gauge theory, although characterized by the most elementary continu- 
ous gauge symmetry, has a rather complex phase structure, induced by the presence 
of monopole excitations. While QED, which shares the same gauge group, appears 
very successfully embedded in the 577(2) x U(l) theory of electroweak interactions, 
it is still quite interesting to study the properties of the pure U(l) lattice gauge 
theory in view of what one can learn from the existence of a non-trivial critical 
point and its implications for possible continuum limits of lattice gauge theories in 
the non-perturbative domain. 

We have devoted several investigations to the study of the properties of the £7(1) 
lattice gauge theory, taking special advantage of a technique we introduced to bridge 
the valley between energy gaps j|, §, and of a new order parameter based on the 
topological properties of monopole networks]!], [jj. Our investigations have been 
based on the study of the model with the Wilson action supplemented by an addi- 
tional monopole term 0] 

S = (3 ]T (l-cos9^) + A^|M p , x | (1.1) 

)1>V,X p,x 

where M P)X = e pa ^ v {Q)^ x+CT - <dp U ,x)/^ and the physical flux e [-it, tx) is 

related to the plaquette angle e (-Arc, Air) by 0^ = 8 W + lim^^ f|. 

It is well known that this system undergoes a phase transition at a coupling f3 CI (\) 
separating a disordered phase, characterized by a condensate of monopole excita- 
tions, for j3 < j3 cr , from an ordered, Coulomb-like phase for j3 > (3 cr . Our own 



2 



studies of the energy distribution have shown that as A increases the phase transi- 
tion changes nature: a gap in the energy distribution indicates that the transition 
is of the first order for A = and for moderately small values of A, whereas the 
disappearance of the gap for larger values of A indicates that the transition becomes 
there of the second order. In this paper we wish to contribute further to the knowl- 
edge of the U(l) lattice gauge theory by exploring in detail the region of larger A 
couplings, where the phase transition appears to be continuous. An important aim 
of our investigation is to verify by a finite size analysis that the transition is indeed 
of second order and to determine deviations of the critical exponent from the Gaus- 
sian case. A further goal is to clarify the behavior of the transition when A becomes 
very large. Because the phase transition is so closely related to the appearance of 
a monopole condensate, exploring in detail the effects of the monopole term in the 
action appears to be of particular interest. 

In this study we use periodic boundary conditions, guaranteeing translational in- 
variance and homogeneity. 

2. The line of phase transitions 

In the paper where they introduced the action ( |1 . 1| ) Barber and Shrock already 
made the observation that there is a shift of the transition point if A is varied [|J. 
It has also been reported in the literature || [7]] that a complete suppression of 
monopoles, corresponding to A = oo, produces the disappearance of the confining 
phase. In our previous investigations we determined the location of the phase tran- 
sition /3 cr (A) for a wide range of values of A Jj], 0. Now we wish to clarify what 
happens at very large A. For this purpose we determine the critical points at val- 
ues of A substantially larger than previously considered and if necessary, allow (3 to 
become negative. 

In order to keep the computational cost for this study within bearable limits we 
have used our topological characterization of the phases [[1], [| to locate the critical 
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points. Our classification is based on the fact that there is an infinite network of 
monopole current lines in the confining phase and no such network in the Coulomb 
phase. On finite lattices "infinite" is to be defined in accordance with the boundary 
conditions ||. For the periodic boundary conditions we are considering here, "infi- 
nite" is equivalent to "topologically nontrivial in all directions" . While for loops the 
topological characterization is straightforward, to determine the existence of an infi- 
nite network it is necessary to perform a more elaborate analysis based on homotopy 
preserving mappings |], 

Because the probability P net to find a network which is nontrivial in all directions 
takes values exactly 1 and in the confining phase and in the Coulomb phase, 
respectively, the frequency of the occurrence of such a network is a very advantageous 
order parameter. Indeed, with an infinite system a single configuration would be 
sufficient to identify the phase. On a finite lattice, because of finite size effects, the 
distinction is not as sharp, but still few configurations are sufficient to discriminate 
bewteen the phases. Examples are given in Figure 1. 

It should be remembered that on finite lattices different order parameters lead to 
slightly different critical (3. On an 8 4 lattice the maximum of the specific heat and 
our topological order parameter give values 1.0075(1) and 1.0074(2) for A = 0, and 
0.3870(5) and 0.372(3) for A = 0.9, respectively. To determine the location of the 
maximum of the specific heat for larger A in an efficient way we first determine the 
critical j3 from the topological order parameter and then find the maximum of the 
specific heat in an easy second step. 

In Figure 2 we show our results for the location of the phase transition /3 cr , defined 
by the maximum of the specific heat, for values of A ranging up to 1.3. It can be 
seen that the phase transition line continues to negative f3. 

Using the topological order parameter we could follow the line of phase transitions 
up to still much larger A: from A = 1.4 where f3 cr = —0.52(2) to A = 10 where f3 cr is 
approximately —1000. It should be emphasized that the characteristic topological 
signatures for the two phases are apparent throughout entire range of A. This shows 
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that, quite remarkably, both phases persist all the way up to very large A. 

We find that finite size effects increase with A. This is indicated by the fact 
that the transition region becomes broader. From Figure 3 it can be seen that the 
width of the peak of the specific heat increases and that its height decreases with 
increasing A. Obviously for very large A a precise determination of the location of 
the maximum of the specific heat becomes quite difficult. The increase of the width 
of the transition region with A is also apparent from the data on the topological 
order parameter P nct in Figure 1 (notice the different scales for f3 in the upper and 
lower parts of the figure). Still, as we already remarked above, the fact that the 
topological order parameter takes values and 1 in the two phases makes it easy to 
discriminate between the two phases even at very large A. 



3. Critical behavior 

In order corroborate our observation that for large A the phase transition becomes 
of second order we have investigated the finite-size scaling behavior of the maximum 
of the specific heat C max . It is expected to be 

Cmax ~ L d (3.2) 

if the phase transition is of first order and 

Cmax ~ L~ (3.3) 

if it is of second order, where a is the critical exponent of the specific heat and v 
the critical exponent of the correlation length. 

In Figure 4 we present results of our simulations for C max , on lattices with L = 6, 
8, 10, 12 for A = 0.9 at the corresponding values of (3 cr . The fit to these data gives 

OL 

- = 0.485(35) (3.4) 
Clearly this is quite far from 4 and thus the transition not of first order. 
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From the hyperscaling relation a = 2 — d v we find 



v = 0.446(5) 



(3.5) 



It is interesting to note that ( |3.5|) is clearly different from the value | of the Gaussian 



case. 



The critical (3 is expected to behave as 



ficr(L) = P cr {oo) + aL 



(3.6) 



From this relation, using the value in Eq. ( |3.5| ) and our data for (3 CT {L) at A = 0.9, 
we get /3 cr (oo) = 0.4059(5) and a = -1.99(6). 

We have also performed simulations for A = 0.8 on lattices with L = 6, 8, 10. It 
turns out that C max does not yet scale in this case. This might be related to the 
particular sensitivity close to a tricritical point. 

In Figure 5 we present the monopole density p = (l/4:L 4 )J2p,x \M PiX \ along the 
transition line. It is interesting to observe that starting at A = 0.9 the density 
becomes constant within the errors of the simulation. 

4. Discussion 

In the past values between 0.33 and 0.50 have been found for the critical exponent 
v, in different contexts |I0]| . This wide range of values has prevented firm con- 
clusions. More precise results, based on a finite size analysis and higher statistics, 
have been recently reported JTTj] for the Wilson action extended by a double charge 
term with coupling 7. 

The well known features of this action are that the first order transition, which 
occurs at 7 = 0, weakens with increasing 7 until it becomes of second order at a 
tricritical point. While for the usual periodic boundary conditions the second order 
region starts at negative 7, the authors of Ref. JTTJ] , using a spherelike lattice, find 
evidence for a second order transition already at 7 = 0. Of course, this could also be 



6 



due to the fact that inhomogeneities weaken the transition M and that the mapping 
from the surface of a hypercube to the 4-sphere can provide inhomogeneities. 

The result v = 0.446(5) presented here and the result v = 0.365(8) given in 
Ref. [JTTJ] are both definitely different from the value \ of the Gaussian case. Thus 
a continuum limit of four-dimensional pure U(l) gauge theory, nontrivial and not 
asymptotically free, appears possible. 

The difference between the two results, though not dramatic, is too large to be 
accounted for by systematic errors. Further studies appear necessary. In particular, 
simulations at some value of A larger than 0.9 will be important for seeing whether 



the value of Eq. (|3.5| ) is universal. Some of the results obtained in Ref. |TT[ for three 
different values of 7 and using different techniques, though close to each other, also 
differ beyond the statistical errors. Results for yet more negative 7 and also for 
periodic boundary conditions would desirable. 

The increase of the finite-size effects with lambda mentioned in Sect. 2 could hint 
at an increase of u, related to a larger correlation length 1] ~ \/3 — f3 cr \~ v '. The 
corresponding decrease of a from 0.216(16) at A = 0.9 to 0, in view of the relation 
C ~ \f3 — /3 CT \~ a , also appears consistent with the behavior of the specific heat seen 
in Figure 3. The behavior of a, in particular, lends support to the fact that the 
increase of finite-size effects is due to the exponents and not to the constants in front 
of the scaling expressions. Possibly v could reach in this way the Gaussian value. 

In regard to our result that both phases are still present for very large A, where f3 cr 
becomes negative, it is to be noted that quite some time ago a similar phenomenon 
has been observed for vortex strings in the 3D 0(2) spin model jL^]. It should 
be further mentioned that negative f3 have also been studied in the compact U(l) 
lattice gauge with complete suppression of monopoles and inclusion of matter fields, 



with the observation of a first order transition at (3 — —0.7 |13|, [IJJ. Clearly this 
is far away from our phase transition line which for very large A runs at extremely 
negative (3. 

We have already shown in Figure 7 of Ref. fl[ that along the phase transition 
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line the average plaquette energy E = (1/6L 4 ) XL >1/a .(l — cos 6^^) increases with 
A. This increase continues for still larger A so that ultimately one comes close to 
the upper limit of the range < E < 2. Thus negative plaquettes contribute 
significantly. In this context it should be noted, however, that, even for A = oo, for 
sufficiently negative (3 the complete suppression of negative plaquettes only reduces 



E by roughly a factor 2 [ 13 ] . 

The non-Gaussian exponent we found implies that the continuum limit of our 
model is different from the usual expectations for QED and would not be seen in 
perturbative approaches. Possibly monopoles survive in the limit. One should also 
be aware of the possibility of obtaining different limits when approaching from the 



strong or the weak coupling regions [|T5 



As a matter of fact, we most likely have a line of critical points in (/?, A) space 
at which different continuum limits may arise if v is not universal. This line is 
to be distinguished from the line of critical points for j3 > (3 cr in case of the pure 
Wilson action, whose existence has been conjectured a long time ago. There the 
scale invariance of Creutz ratios fl6j hints at a diverging correlation length and 
thus at a line of critical points. Similar indications come from calculations of the 
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Further work appears necessary to fully explore the critical properties of the model 
and, in particular, to determine whether its critical points can give origin to a viable 
theory. A related question is which form of the action is the most profitable to 
use. The addition of a monopole term is attractive because of the close relation 
of monopoles to the phase structure. We wish to note, however, that instead of 
^ J2 px \M PtX \, nJ2 px (.M PtX ) 2 or some other form might turn out to give better results. 
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Figure captions 



Fig. 1. Order parameter P net as function of (3 for A = and A = 0.9 
on an 8 4 lattice. 

Fig. 2. Location of the phase transition point (3 CT as function of A 

for 8 4 (squares) and 16 4 (crosses) lattices. 
Fig. 3. Specific heat C as function of (5 for A = 0, 0.9, 1.3 on an 8 4 lattice. 
Fig. 4. C max versus L for A = 0.9 at f3 cr . 
Fig. 5. Monopole density as function of A. 
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